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The conventional approach to the turbulent energy cascade, based on Richardson-Kolmogorov 
phenomenology, ignores the topology of emerging vortices, which is related to the helicity of the 
turbulent flow. It is generally believed that helicity can play a significant role in turbulent systems, 
e.g., supporting the generation of large-scale magnetic fields, but its impact on the energy cascade 
to small scales has never been observed. We suggest for the first time a generalized phenomenology 
for isotropic turbulence with an arbitrary spectral distribution of the helicity. We discuss various 
scenarios of direct turbulent cascades with new helicity effect, which can be interpreted as a hindering 
of the spectral energy transfer. Therefore the energy is accumulated and redistributed so that the 
efficiency of non-linear interactions will be sufficient to provide a constant energy flux. We confirm 
our phenomenology by high Reynolds number numerical simulations based on a shell model of helical 
turbulence. The energy in our model is injected at a certain large scale only, whereas the source of 
helicity is distributed over all scales. In particular, we found that the helical bottleneck effect can 
appear in the inertial interval of the energy spectrum. 

PACS numbers: 47.27.Gs, 47.27.E- 


Many years ago, K. Moffat delivered a verdict on the 
influence of helicity on spectral properties of turbulent 
flows Q: 

no matter how strong the level of helicity in¬ 
jection may be at wave-numbers of order kg, 
the relative level of helicity as measured by the 
dimensionless ratio H(k)/2kE{k) must grow 
progressively weaker with increasing k; and 
when k/l^ is sufficiently large it may be con¬ 
jectured jjy that the helicity has negligible dy¬ 
namical effect, and is itself connected and dif¬ 
fused in much the same way as a dynamically 
passive scalar contaminant 0/. 

Thus the mean helicity can play a significant role in tur¬ 
bulent systems, e.g., supporting the generation of large- 
scale magnetic fields [1, @ , but its impact on the spectral 
properties of turbulent flows is practically negligible. 

The experimental study of helicity needs hardly afford¬ 
able three-dimensional measurements with high spacial 
resolution. Moreover, the helical forcing of a turbulent 
flow is not straightforward to implement in the labora¬ 
tory. This is why the impact of helicity on the energy 
cascade has been studied in few laboratory experiments. 
The energy and helicity spectra with a slope ”-5/3” have 
been measured in the atmospheric boundary layer Q. 

Direct numerical simulations are able to model helical 
isotropic turbulence and evaluate the spectral properties 
of the flow. They confirm the existence of Kolmogorov 
spectra (see for example 0, S|)- The same results have 
been demonstrated at a high spectral resolution with the 
help of the shell model of turbulence 00, showing no 
reason to doubt Moffat’s argument. 

Features of helical turbulence are considered able to 
address the problem of the inverse energy cascade. Like 


two-dimensional turbulence, it was suggested that two 
inertial ranges could be realized corresponding to the di¬ 
rect cascade of helicity with an energy spectrum ”-7/3” 
and an inverse cascade of the energy with a slope ”-5/3” 
iEl. However, this analogy is not simply applicable 
because the helicity spectral density, in contrast to en- 
strophy, is not proportional to the energy spectral den¬ 
sity. Recently, direct numerical simulations of the three- 
dimensional isotropic turbulence confirmed that the in¬ 
verse cascade appears in truncated three-dimensional 
Navier-Stokes equations 12|, in which the helicity is pre¬ 
scribed to be a sign-defined quantity and therefore di¬ 
rectly related to the energy. 

However, the statistical properties of helical turbulence 
are an issue. Indeed in some cases the_spectra of turbu¬ 
lent magnetohydrodynamical (see [ISEl and references 
therein), convective [l5j and atmospheric flows [l^ differ 
from Kolmogorov’s scaling, which can be interpreted as 
an effect of helicity. Furthermore, the influence of helicity 
on the flow properties can be manifested in the behavior 
of the turbulent spectrum of a passive scalar (e.g. an 
aerosol) 0. Helicity can change not only the energy 
and passive scalar spectra, but can lead to the reduc¬ 
tion of the turbulent viscosity and energy accumulation 
at large scales 18[, or can play a significant role in the 
generation of large-scale helical vortex structures in the 
atmospheres of the Earth and planets (l9l - [^ . To clarify 
the role of helicity, we consider the behavior of helical 
isotropic turbulence in a somewhat different way. 

Coming back to forced helical turbulence, we revise 
the view that helicity is injected into the flow together 
with energy at the same scale. Theoretically, one can 
assume that turbulence is excited by a source of energy 
at the largest scale and an independent source of pure 
helicity, acting at a certain scale or over all scales in the 
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inertial interval. Then, the helicity spectral flux is not 
constant anymore and the helicity spectral density can 
reach significantly higher values and influence the energy 
cascade. Real physical situations usually are far from 
ideal, but can be similar to some extent e.g. in rotating 
convective flows [ 2 ^ . 

In the present work, we aim to elaborate on a phe¬ 
nomenology of an energy spectrum for the highly helical 
turbulence. Our theoretical result is justified by numer¬ 
ical simulations of helical turbulence with a controlled 
constant energy flux and distributed helicity injection 
over the scales in the inertial range. We demonstrate 
various scenarios of turbulent cascades with the helicity 
effect, and discuss its physical meaning with regard to 
the bottleneck phenomenon known for conventional tur¬ 
bulence 


24|. 


First, we adopt the basic statement of Kolmogorov’s 
approach, which claims that in the inertial range at any 
scale n the energy flux is equal to the dissipation rate, 
= £. We consider a geometric sequence of scales and 
corresponding wave numbers l~^ ^ kn ^ A". The energy 
flux is related to the velocity pulsations at this scale 
as 

nf R:! vlkn (1) 

and the energy of velocity pulsations En ^ ^ 

{e/knY^^- 

Second, we follow the decomposition of velocity pul¬ 
sations in two helical modes [llj, -I- with 

corresponding energies ^ Then the energy 

and helicity at the scale n are 


En — E^ + E^ , 
H^ = H+ + H-= k^{E+ - E-). 


( 2 ) 


The energy flux at scale n is decomposed into four 


terms: 


kn{v+f and kn{v^f 


The contribution of the first two terms dominates in local 
interactions in the inertial range and becomes compara¬ 
ble with the two other terms for non-local interaction 
[111. The significant role of terms kniv^Y distin¬ 
guished in the inverse cascade [l^. For the direct cas¬ 
cades of the energy and helicity one can neglect the con¬ 
tribution of these terms. 

We introduce the relative helicity = Hn/{knEn), 
which allows us to link the intensity of two helical modes 


- + n-H- 

" i + m 


~ ■ 


Then the energy flux provided by local interactions can 
be estimated as 


Replacing (HD by © we finally obtain 

2/3 

Ef ' ' 


( 4 ) 


kni^n + Cn) 


( 5 ) 


One can express E^ = ^^E^ from i© and obtain the 
total energy 


E„ = EY + E- (eCn/knf'' : 
where the dimensionless variable 


Cn = 




( 6 ) 


( 7 ) 


depends on Hr{k) and defines “the degree of helical 
blocking” of the spectral energy flux at a given scale. 
\Hr{k)\ characterizes the dominance of some helical 
modes over others with the opposite sign, i.e. it is 
the helical part of the energy. Then a new parameter 
S{k) = 1 — \Hr{k)\ corresponds to the non-helical part of 
the energy, which is free of helicity. For the highly helical 
case iF” ±1, formula (13 has as asymptote 


The corresponding spectral energy density 

E{k) « £2/3A:-5/35(fc)-i/3 


( 8 ) 


( 9 ) 


is independent of the sign of the injected helicity. Usually 
for the single-scale forcing of helicity, Ff” oc k~^ and the 
parameter does not differ from unity. One can expect 
a significant change in the turbulent spectra for highly 
helical turbulence only. 

We interpret our results as a hindering effect of the 
spectral energy transfer in scales with high relative helic¬ 
ity. The energy should be accumulated and redistributed 
so that the efficiency of non-linear interactions will be 
enough to provide the constant energy flux, which is pre¬ 
determined by the energy injection rate. A similar con¬ 
sequence is observed as a result of the bottleneck phe¬ 
nomenon 2M the non-helical turbulent cascade when 


non-local interactions drop out of the spectral energy 
transfer at the end of the inertial range. We exploit this 
analogy to name our effect the helical bottleneck effect, 
having in mind the helical mechanism of cascade block¬ 
ing. 

Further study of the suggested helical effect W DNS 
would be highly desirable. Recent simulations 2^ con¬ 
firm a change of energy spectrum due to distributed 
helicity and justify a dominant contribution of triads 
(3) kn{vYYvi^ which is assumed in Eq. ©. However ver¬ 
ification of Eq. (© is hardly affordable by DNS with a 
resolution which seems to be achievable in near future. 
We demonstrate the realizability of the energy spectrum 
(© with helical correction dzD by numerical simulation us¬ 
ing the helical shell model of turbulence (HI type model 
as it is classified in the review ). The governing equa¬ 
tions for complex shell variables C/„ are 


iln = Wn{U, U) - Re-^klUn + /,f + /, 


H 
n ’ 


( 10 ) 
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FIG. 1: Energy {E(k), points) and normalized Felicity {H{k)/k, circles) spectra compensated by (Kolmogorov’s slope 

with intermittency correction) are shown in a grid of panels for each pairs of parameters 770 = 0.3, 1, 3.32 (from top to bottom) 
and a = 0,0.5, 1 (from left to right). Exceptional cases are panel (c) 770 = 0.03, a = 1.75 and panel (g) 770 = 10, a = —0.5. 
Theoretically predicted helical correction to the energy spectrum (((fc) (defined by Eq. ((7|) is shown by the solid line. 


where = A" is the wave number attributed to the shell 
n, A characterizes the shell thickness, a bilinear operator 
Wn determines interactions of with its neighbors Un+i 
and Un-i- The energy forcing term is prescribed to 
operate at the largest scale only, i.e. in the shell n = 0 
{ko = 1), and provides the mean energy injection rate e 
and the zero mean helicity injection rate. For the helicity 
forcing, we introduce a scale distributed force 

/f = imKUniul + c/f )/2, (11) 

which acts in all shells and includes two parameters r/o 
and a (the asterisk is for complex conjugation), r/o 
defines the helicity injection rate at the largest scale 
(tt = 0 ), which is the energy forcing scale also, a specifies 
the scaling of the helicity influx. 

Variables C/„ define the energy En = |t/„p/2 and helic¬ 
ity Hn = ikn{U*^ — of the shell n. The correspond¬ 
ing spectra are E{kn) = En/kn and H{kn) = As 

decomposition m, one can split En and Hn into two he¬ 
lical parts. Then the helicity injection rate at the shell 
n is Hn = 4?7ofc““''^A+i?“. Note, that the force fn satu¬ 
rates and tends to zero if the helicity Hn approaches the 


limiting value fc„|C/jip. 

Numerical simulations of shell model equations m 
are performed for Re = 10® and A = 1.62. Statistically 
stationary spectra and fluxes are obtained by averaging 
over 256 runs (with random initial conditions and en¬ 
ergy forcing), lO’^ time units long each. We choose the 
combination of governing parameters 770 = 0.3,1,3.32 
and a = 0,0.5,1 to study the transition of turbulence 
from low to high intensity of helicity injection (controlled 
by rjo) and from mostly large-scale injection to injection 
equally distributed over all scales (controlled by a). Two 
special cases with r/o = 0.03, a = 1.75 and ryo = 10, 
a = —0.5 are included to emphasize the effect of he¬ 
licity. The corresponding energy E(k) and normalized 
helicity H{k)/k spectra are presented in Fig. [TJ Both 
spectra are compensated by Kolmogorov’s slope ”-5/3” 
to highlight the difference in comparison with conven¬ 
tional turbulence. To get the flat energy and helicity 
spectra, we divide by a slightly steeper slope (caused by 
intermittency), namely, by One can compare the 

result of simulations with the theoretically predicted en¬ 
ergy spectrum which is shown by a solid line for each 
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set of the helical forcing parameters. 

Panels (a) and (b) demonstrate the cases of weak he- 
licity injection (779 = 0.3) in which passive helicity does 
not change energy spectrum. The helicity spectrum at 
q; = 0 scales as k~^, which looks like the spectrum for a 
passive scalar. The moderate increase of a increases the 
slope of the helicity spectrum only (panel b). 

Panels (d), (e) and (f) show results for stronger helic¬ 
ity injection (770 = 1) at which the helicity is becoming 
active. Depending on a, the energy accumulates in cer¬ 
tain shells at large scales, i.e. in only four shells (panel d) 
or in all shells with < 100 (panel f). In these ranges, 
the relative helicity tends to unity, which can be dis¬ 
tinguished since H{k)/k (open circles) approaches E{k) 
(closed circles). 

Panels (h) and (i) correspond to a further gradual in¬ 
tensification of helicity injection (770 = 3.32), which leads 
to a steeper slope of the energy spectrum. Again de¬ 
pending on a, this can be observed in a part (panel h) 
or almost the entire inertial range (panel i). In the latter 
case, we find E{k) ~ over three decades of scales. 

Panel (g) presents a special case with the largest 770 = 
10 and lowest a = —0.5. Two inertial ranges can be 
recognized with pronounced active and passive behavior 
of helicity. A completely helical turbulent cascade, that 
is E[k) - and H[k) - appears at k < 10 . 

Kolmogorov’s slope is recovered in the rest of the scales 
and there is a sharp border between the two ranges. 

Panel (c) shows the opposite case with the lowest 
770 = 0.03 and largest a = 1.75. This combination of the 
forcing parameters provides a passive cascade of the helic¬ 
ity at the large scale and an active one at the small scale. 
This active range of scales is characterized by the high 
relative helicity, which decelerates the energy cascade and 
results in the pileup of energy near the small-scale limit 
of the inertial range. This regime exaggerates the realis¬ 
tic scenario of small-scale helicity injection. However, it 
has particular interest because the result is similar to the 
bottleneck effect 2J, [27| . Note, that spectra calculated 
using theoretical formula 0 , which are shown by solid 
lines in Fig. [T] fit well the numerical spectra in all cases. 

The basic feature of the inertial range dynamics is the 
constant spectral energy flux. Fig. Ha) presents the spec¬ 
tral fluxes of energy and helicity for three cases. The 
energy fluxes are not disturbed by the forcing m , even 
for a very strong injection of helicity. Apparently the 
helicity spectral flux is influenced by distributed helic¬ 
ity forcing. The helicity spectral flux continually ap¬ 
proaches the asymptote k with an increase of 770 and a. 
n^(fc) - kU^{k) is the theoretical limit, which corre¬ 
sponds to the maximally helical turbulent cascade. 

Instead of the relative helicity Hr{k), the final expres¬ 
sion ([9]) includes S{k). Both quantities Hr{k) and d{k) 
indicate scales where the energy is affected by the helic¬ 
ity {Hr{k) > 0.7 or d{k) < 0.3). But S{k) has a constant 
slope if the energy spectrum obeys a power law in some 




FIG. 2: (a) Spectral fluxes of energy and helicity (inset), (b) 
Distributions of S{k) = 1 — \Hr{k)\. 


range of scales (see Fig. Hb))- Therefore, we conclude 
that the ’’free” part of the energy 5{k) is the more ap¬ 
propriate quantity which controls the efficiency of the 
spectral energy transfer. 

Following Eq. we should expect S{k) ~ if the 
power law "-1 jZ" appears in the energy spectrum. This 
is too steep to be realized over the whole inertial range. 
The limited range of scales with one order of magnitude 
is more realistic, as we find for the case in Fig.[T](g)- The 
affected range of scales shortens as the helical correction 
strengthens. This is another argument in favour of the 
analogy with the bottleneck effect. 

Our results expand the common point of view that the 
helicity is a possible trigger for inverse energy transfer, 
which support long-living large-scale vortical structures 
in astrophysical and geophysical flows [ 2 ^. We find a 
more modest role for helicity in reducing the efficiency 
of spectral energy transfer to small scales. There is no 
inverse cascade, just energy accumulation due to the hin¬ 
dered cascade in the highly helical flow. For example, we 
suggest that the helical bottleneck effect could be respon¬ 
sible for steepening the experimentally obtained energy 
spectrum in a von Karman swirling flow 23|. In that 
experiment, the helicity was intensively injected by the 
blades at smaller scales than the scale of energy injec¬ 
tion by rotating discs [ 2 ^. Another applicable situation 
might be a flow with shear generated helicity, which is 
favorable for development and intensification of tropical 
cyclones [30. We note that the spacial imbalance 
of helicity can support a local intensification of vortical 
structures, even if the global helicity is negligible. 

We revise the conventional approach to the turbulent 
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energy cascade, based on Richardson-Kolmogorov phe¬ 
nomenology, where the topology of emerging vortices is 
ignored. This new generalization of the Kolmogorov phe¬ 
nomenology takes into account the arbitrary spectral dis¬ 
tribution of helicity and describes a deviation from the 
-5/3 power law for the energy spectrum. The immediate 
consequence for the mean field theory is that an esti¬ 
mate of the effective turbulent diffusivity has to consider 
the helicity. This is an important message for physicists 
dealing with sub-grid models of turbulence. 
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Academy of Science, Ural Branch. 
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